Importance of carbon solubility and wetting properties of nickel nanoparticles for 

single wall nanotube growth 
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Optimized growth of Single Wall Carbon Nanotubes requires a full knowledge of the actual state 
of the catalyst nanoparticle and its interface with the tube. Using Tight Binding based atomistic 
computer simulations, we calculate carbon adsorption isotherms on nanoparticles of nickel, a typical 
catalyst, and show that carbon solubility increases for smaller nanoparticles that are either molten or 
surface molten under experimental conditions. Increasing carbon content favors the dewetting of Ni 
nanoparticles with respect to sp 2 carbon walls, a necessary property to limit catalyst encapsulation 
and deactivation. Grand Canonical Monte Carlo simulations of the growth of tube embryos show 
that wetting properties of the nanoparticles, controlled by carbon solubility, are of fundamental 
importance to enable the growth, shedding a new light on the growth mechanisms. 

PACS numbers: 61.46.Fg; 68.37.Lp 
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The currently most widespread synthesis technique 
for single wall carbon nanotubes (SWNT) is the Cat- 
alytic Chemical Vapor Deposition (CCVD). Although 
massively developed in the last decade, the process still 
gives rise to unsolved questions that hamper the much 
coveted progress towards a selective growth of SWNTs, 
with designed chirality and, hence, properties. We fo- 
cus here on two questions of importance in the field of 
nanoscience that have far reaching consequences for nan- 
otube growth. In CCVD, for SWNT as well as graphene 
synthesis, a carbon rich precursor interacts with a metal- 
lic catalyst surface, possibly leading to carbon dissolution 
and structural modifications. 

How carbon solubility and physical state of the 
nanoparticle (NP) depend on temperature and size is 
the first question we address. It has been empirically 
shown that efficient catalysts for SWNT growth display 
a limited, but non zero carbon solubility [lj]. Ding et 
al. could calculate the melting temperatures of Fe- 
C NPs with given size and composition. Under the as- 
sumption that smaller NP size would increase the inter- 
nal (Laplace) pressure, Harutyunyan et al. 0] concluded 
that smaller Fe-C NPs display a lower carbon solubility 
limit. For other transition metal-carbon systems (Ni-C, 
Pd-C), subsurface interstitial sites have been identified 
as most favorable for carbon incorporation jH-Q] on flat 
surfaces. Assuming that this still holds true for NPs, one 
could expect the opposite result that carbon solubility in- 
creases for smaller NPs of these metals, because a large 
fraction of C should be expected close to the surface and 
the surface to volume ratio becomes larger for smaller 
sizes. In order to get a direct answer, we calculate the 
carbon adsorption isotherms on nickel NPs. 

The second aspect is to understand how this modified 
carbon solubility influences the wetting and interfacial 
properties of nickel NPs in contact with carbon sp 2 lay- 



ers. Metals are generally reported not to wet carbon 
nanotubes but capillary effects can favor the pen- 
etration of metal catalyst inside the tube In their 
study of the catalyst restructuring during SWNT growth, 
Moseler et al. [lfj assumed a 180 ° contact angle of Ni 
with graphene. As also done by Schebarchov et al. [§] 



and Borjesson et al. 11|, the nanotube wall was consid- 



ered rigid in their Molecular Dynamics simulation, and 
no carbon solubility was taken into account. However, 
back in 1971, Naidich et al. [H] showed that the con- 
tact angle 9 of a macroscopic Ni (and also Co and Fe) 
drop on graphite is strongly modified by C dissolution. If 
this strong change of wetting properties of these typical 
catalysts holds true at the nanoscalc, it certainly plays 
a major role in the growth mechanisms of SWNTs, e.g., 
to prevent the encapsulation of the NP by the growing 
tube wall. To study this, we directly simulate the struc- 
ture of carbon enriched nickel NPs on a graphene layer 
at different compositions and temperatures. In a last 
part, we show that the enhanced carbon solubility and 
dewetting properties of the NPs critically influence the 
growth of the tube wall, that is reported to have a major 
effect on the tube chiral distribution [H, 14|, and there- 
fore evidence an essential, though somewhat overlooked, 
property of the metal catalyst. 

To perform these calculations, we use a carefully val- 
idated tight binding model to describe the interactions 
between nickel and carbon atoms 15|. This model is in- 
cluded in a Monte Carlo (MC) code, working in the grand 
canonical (GC) ensemble, which is the natural scheme to 
study the incorporation of carbon in a simulation box, 
with given chemical potential (/ic) and temperature (T). 
The technical implementation of the model has recently 
been made very efficient enabling us to extend the 
range and accuracy of curious calculations 0. To 
simulate CCVD, that is a surface reaction, C atoms are 
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FIG. 1: Structures of C-rich NPs (405 Ni; C fraction cor- 
responding to the solubility limit; cut through the NP; Ni 
atoms: orange, C atoms: grey), and analysis using the distri- 
bution of radial atomic densities (black : C, red : Ni) from 
a central Ni atom (the closest from the NP center of mass). 
Bottom panel: 1400 K, particle is molten, C atoms are dis- 
tributed in the whole NP. Top panel: 800 K, the core-shell 
structure is clearly visible, with a molten outer layer with C 
atoms, dimers and some trimers, while the crystalline core is 
pure Ni. 



tentatively added or removed in a "peel" of about 0.25 
nm thickness from the NP surface. As explained in [l8| 
this algorithm drives the system towards its equilibrium 
at fixed N/vi, [ic an d T. Our model yields a melting 
temperature for pure bulk Ni (T m (bulk) = 2010 ± 35 K) 
that is ~15% too high, while the relative melting tem- 
peratures of pure NPs (T m (NP) / T m (bulk) are 0.84, 
0.81, 0.76, 0.69, 0.59 for NPs with 807, 405, 201, 147 and 
55 Ni atoms, respectively [l9j. We could also check that 
the maximum carbon solubility in a bulk system of 576 
Ni atoms never exceeds 4-5% in the crystalline phase, in 
agreement with the experimental phase diagram. 

How carbon solubility is modified by the nanometric 
size of the Ni particles can be studied by computing the 
carbon adsorption (or incorporation) isotherms on NPs of 
various sizes. We start from a bare Ni NP and run a series 
of GCMC simulations at different \ic and T. Computa- 
tional details are given in Supplemental Material. Once 
equilibrium is reached, we record the average number of 
C adsorbed inside and possibly outside the NP. Carbon 
inside the NP generally appears as dimers (C-C distance 
around 0.19 nm) close to the surface. This is exemplified 
in Fig. [T] that presents two NPs with 405 Ni, saturated 
with carbon, at 800 K and 1400 K. The former displays 
a crystalline core and a disordered outer shell with some 
C dissolved, while the latter is molten. Calculating the 
"molten" fraction is readily done using pair correlation 
functions. We see that the crystalline core contains few 
carbons, while the external molten shell is C-saturated in 
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FIG. 2: Carbon adsorption isotherms calculated at 1000 K 
for nickel particle sizes 55, 147, 201, 405 and 807 atoms. 55 
and 147 particles are icosahedral when pure, the others have 
the equilibrium Wulff shape of the finite face centered cubic 
structure. The y-axis is the carbon content inside the parti- 
cle. Dotted lines indicate that, at the next \ic step, carbon 
atoms appear on the surface. The inset presents the relative 
thickness of the outer liquid layer (p), as a function of \xc- 
NPs with 55 and 147 Ni atoms are liquid as soon as carbon is 
adsorbed, in these temperature and carbon chemical potential 
conditions. 

these examples. "Molten" here means disordered, since, 
using a MC method, we have no time scale and hence 
no dynamics and diffusion coefficients: we cannot distin- 
guish liquid and amorphous structures. Fig. [2] presents 
the average mole fraction of carbon dissolved inside NPs 
as a function of \ic at 1000 K. Data are plotted up to 
the carbon solubility limit only, defined as the carbon 
fraction beyond which carbon species (dimers, chains ...) 
begin to appear on the surface of the NP. The general 
trend is that carbon incorporation begins at lower /ic for 
smaller NPs. At 800 and 1000 K though, an additional 
effect of the NP structure can be identified. NPs with 55 
and 147 are icosahedral and present a close packed sur- 
face that makes carbon adsorption or incorporation less 
favorable than on the (100) facets or edges of the Wulff 
shaped particles with 201, 405 and 807 Ni atoms. This 
explains why the adsorption curve for a M147 icosahedral 
NP plotted in Fig. [2J lies below that for a Ni 20 i Wulff 
shaped one, as long as the NPs are not molten. Another 
way to read this graph is to notice that at a given carbon 
chemical potential, smaller particles display larger car- 
bon content. The incorporation of carbon modifies the 
physical state of the NP. At 1000 K, all particles with 
55 and 147 Ni atoms are molten or amorphous, as soon 
as C is incorporated. Depending on nc, and hence on 
their carbon content, larger particles (201 to 807 Ni) can 
be either molten, or display a core shell structure. The 
crystalline fraction of the NPs can be deduced from the 
inset of Fig. [3J and a full phase diagram can in prin- 
ciple be computed for NPs of different sizes. Different 
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from the conclusions of Harutyunyan et al. 0] obtained 
for the Fe-C system, we show here that in the case of 
Ni-C NPs, the solubility is larger when the particle size 
is smaller. This difference stems from the fact that calcu- 
lations reported in Ref. Q rely on the assumption that 
smaller size induces a larger Laplace pressure within the 
NP. This is true for a pure NP, but our calculations show 
that adsorbed C atoms strongly interact with the Ni sur- 
face and modify the surface properties of the NP, and 
hence its internal pressure. 

In the size range (1-3 nm) studied here, that is typi- 
cal for SWNTs growth, the NP surface is strongly dis- 
ordered by C incorporation and fairly smooth. TEM 
observations, performed in situ on larger NPs HH 
or after growth [2^, often show carbon walls ending on 
step edges of the metal catalyst. These step edges are 
sometimes considered as necessary for growth. Our cal- 
culations show that growth can take place without any 
step edge and that the movement of the tube with re- 
spect to the NP results from another driving force. In 
situ Transmission Electron Microscopy (TEM) images of 
carbon nanofibers growth [20I [23| , as well as theoretical 
investigations 0, [l(j] strongly suggest that the dewetting 
tendency of the NP and the tube should play a role. We 
therefore study the wetting of such NPs in contact with 
sp 2 carbon layers. 

To separate the capillary absorption that can occur 
at these small nanotube sizes Q from pure wetting, we 
calculate (Fig. [3]) the contact angle of NPs with 405 Ni 
atoms and C fractions from zero to the solubility limit, 
close to 25%, on a graphene layer of 3840 C atoms, us- 
ing canonical Monte Carlo simulations. A visual inspec- 
tion of the resulting structures (Fig. [3|) reveals a qualita- 
tive agreement with the experimental data [12[ obtained 
for macroscopic droplet sizes: the contact angle becomes 
larger and beyond 90 ° when carbon fraction increases. 
The accuracy of our calculations is limited to about ± 
15 ° . This results more from the non-symmetric droplets 
shapes (caused by their very small sizes and the atomic 
roughness of the substrate), than from statistical errors 
over the hundred configurations considered for averag- 
ing. However, the trend is clear and, either because of 
the nanometric size of the NP, or because of the interac- 
tion with one single graphene layer, or both, the contact 
angles are larger than reported in Naidich et al. [IH for 
a macroscopic sample. 

In order to see how the SWNT growth mechanisms are 
influenced by this modification of the wetting properties 
of the NP when C is dissolved in, we devised the fol- 
lowing set of computer experiments. We start by fixing 
short nanotube butts with about 80 C atoms on pure Ni 
NPs with 70 or 85 atoms. The tubes are chosen with 
diameters in the 0.7-0.8 nm range and different chiral- 
ities. These atomic configurations are first moderately 
heated and subsequently relaxed to make them reach 
their state of minimal energy: some Ni atoms tend to 
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FIG. 3: Contact angles of a N1405 NP on graphene, as a func- 
tion of the carbon fraction, at 1000 (red) and 1400 K (green) . 
The two snapshots correspond to pure (left) or carbon satu- 
rated (right) NPs. 



enter inside the capped tube, while remaining Ni atoms 
stay outside. We used them as starting configurations 
for a series of GCMC simulations with temperature and 
fic conditions chosen to allow the growth of the tubes. 
Out of more than a hundred trials spanning a large pa- 
rameters range, we selected 10 runs for which growth 
was effectively promoted, for around -6.5 cV/atom 
and 1050 K <T<1300 K. For all samples, the growth 
starts only once the carbon content of the NP has reached 
about 15%. This is slightly below the equilibrium car- 
bon concentration values calculated above for bare NPs 
(around 20%), because the C distribution within the NPs 
is slightly inhomogeneous, lower in the part of the NP 
that is protected by the carbon cap that prevents C ad- 
sorption. Attempted carbon insertions are accepted only 
on the free surface of the NP and, once it is saturated, 
C atoms diffuse on the surface to form very mobile short 
chains that eventually get attached to the tube lip like a 
garland. Growth by attachment of short carbon chains 
has been evidenced for graphene on metal p4| , and chain 
length probably depends on the experimental conditions. 
Fig. |4] presents a series of atomic configurations during 
a successful growth. As outlined in our recent combined 
TEM and computer simulation study [25| , discriminating 
so-called perpendicular and tangential growth modes, the 
most efficient one is the latter: it corresponds to the sit- 
uation depicted here, with carbon walls developing along 
the NP's side and not pushing perpendicular to the NP 
surface. The growth is measured by the distance from the 
topmost C atom of the initial butt to its closest Ni atom 
and by the number of carbon rings formed, hexagons, but 
also defects such as pentagons or heptagons. In addition, 
we see the sidewall growing along the C-enriched Ni NP 
up to MC cycle 300. We then notice a very fast extru- 
sion of the NP from the growing tube. Starting from a 
pure Ni NP, all successful growths proceed in this way. 
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FIG. 4: Top panel : Snapshots of the atomic configurations 
during the growth of a tube butt (80 C atoms: grey) on a NP 
(85 Ni atoms: yellow), a) Initial configuration as prepared, b) 
after moderate heating and subsequent relaxation, this con- 
figuration is used to start GCMC calculations, c) MC cycle 
300, before cap detachment, added C atoms are in blue, d) 
MC cycle 350, after cap detachment, d) MC cycle 450, during 
growth. Lower panel, left: growth length d (black, defined 
above) and mole fraction of C dissolved in the NP (red), as a 
function of the Monte Carlo cycles. Right: number of carbon 
rings formed (blue: hexagons, black: pentagons, green: hep- 
tagons) and number of added C atoms (red), as a function of 
the Monte Carlo cycles. 



The delay between the moment when the critical carbon 
concentration is reached and the actual cap detachment 
is reminiscent of the contact angle hysteresis observed for 
growing or receding droplets on a surface: roughness, the 
amount of NP surface covered by sp 2 carbon, as well as 
the interaction between the tube lip and the NP, certainly 
plays a role. It also frequently happens, at low tempera- 
tures and high fic favoring C incorporation, that the cap 
does not detach and the growth of the side wall results 
in a particle encapsulation. This dewetting tendency is 
reversible: if we remove the dissolved C atoms, the NP 
tends to go back inside the tube (see Supplemental Ma- 
terial). 

This analysis enables us to suggest that the SWNT 
growth results from two competing phenomena. On the 
one hand, once a carbon cap is nucleated on the NP that 
serves as a template, the amount of C dissolved in the 
NP favors dewetting conditions that tend to separate the 
growing tube from the NP. On the other hand, the carbon 
wall growing alongside the NP and the interaction of the 
NP with the tube lip prevent a complete detachment. In 
steady state regime, these two moves are balanced and a 
continuous growth is observed, with a global kinetics gov- 
erned by the corresponding kinetics of these two moves, 
plus the feedstock decomposition and the carbon diffu- 
sion kinetics. The last two processes do not depend on 
the tube chirality. Both remaining processes, tube dewet- 
ting from the NP and wall growth through carbon chains 
attachment at the tube lip strongly depend on the tube 



- NP interaction and are probably a key to chiral se- 
lectivity, giving a clue to interpret recent experimental 
findings. We note that Motta et al. [26| report improved 
growth, yielding longer tubes, when sulfur is included in 
the feedstock, acting as a surfactant. Sundaram et al. [27} 
reported the synthesis of metallic SWNT fibers suggest- 
ing that sulfur limits the growth of the catalyst NPs. 
Harutyunyan et al. [2S| observed a preferential growth 
of metallic SWNTs by modifying the ambient gas phase, 
and explained it by a modification of the NPs morphol- 
ogy induced by adsorbed species, while Chiang et al. [29[ 
could link the chirality distribution to the composition of 
their NixFe^-^) catalyst. From our analysis, we suggest 
that one could also interpret these results as a modifica- 
tion of the interfacial properties of the metal NPs, that 
depend very much on their size, the temperature and the 
carbon chemical potential that control the carbon con- 
tent. 

To conclude, our findings can be summarized as fol- 
lows. We first showed that in the case of nickel that 
presents a moderate phase separation tendency with car- 
bon, carbon solubility tends to increase when the NP 
size becomes smaller. Carbon atoms are concentrated 
close to the surface, often as dimers. Depending on their 
size, carbon content and temperature, NPs can be either 
molten or partially molten with a crystalline core. We 
then showed that the increase of the contact angle with 
increasing C content, already noticed for a macroscopic 
Ni drop on graphite, is even more pronounced for a Ni NP 
on graphene, with contact angles varying from 80 ° (pure 
Ni) to almost 170 ° for C-saturated Ni. Since this C in- 
duced dewetting property has been reported for Ni, Co 
and Fe that are typical catalyst for SWNT growth, it 
is very likely that it should be included in the portfo- 
lio of necessary properties for a good catalyst. Lastly, 
using GCMC calculations starting from a preformed car- 
bon tube butt fixed on a pure Ni NP, we showed that 
the SWNT growth takes place preferably by carbon wall 
growing along the NP side, the dissolved C promoting 
the dewetting of the NP from the tube, and also causing 
the NP surface to be strongly disordered (liquid or amor- 
phous) . We suggest that engineering the reactivity of the 
tube lip in the presence of a disordered metal, enriched 
with carbon, and the interfacial properties of the NP are 
keys to chiral selectivity. 
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